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Abstract 

Comonotonicity had been a extreme case of dependency between ran- 
dom variables. This article consider an extension of single life model 
under multiple dependent decrement causes to the case of comonotonic 
group-life. 

1 Introduction 

The concept of comonotonicity was proposed for research of economics, and it 
had been introduced into the insurance actuarial in 1990's (see. Dhaene J. et 
aL[3]) for ordering risks and stop- loss premium determination (see. Dhaene J. 
et al. [3]). Research results show that comonotonicity had been a extreme case 
of dependency between random variables (see. Dhaene J. et al. [3] , [4] , [5] , Ribas 
C.[5]). Thus, for example, comonotonicity between lives in a group implies a 
special characteristics of the group-life that one member of the group attains 
venerable age then all of other members of the group also attain venerable age. 
The decrement models with respect to multiple dependent causes of withdrawal 
have been received a attention from life insurance area, Carriere JF. (1994, 1998) 
[l],[2] and Vladimir K. et aZ.(2007)[9] consider the survival of the single life 
under multiple dependent decrement causes, their results are mainly based on 
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the copula functions. A statistical estimation on the net survival functions us- 
ing the observations of crude survivals was proposed. Wang et aZ.(2009) [TO] 
has considered the decrement model for comonotonic group-life with status of 
joint-survival and last survival under independent multiple decrement causes. 
The distribution model with statistical analysis on comonotonic group- life under 
multiple dependent decrement causes will be considered in the present paper. 
The comonotonicity properties seemingly make that each group-life can be rep- 
resented by any single member of the group, which would leads to a simple way 
inducing a joint distribution for the considered collection of comonotonic group- 
life. This paper is organized as follows. In Section 2, the main concepts like 
comonotonicity, copula functions, survival function as well as some useful ob- 
tained results are mentioned, in Section 3, we construct a model for dependent 
multiple comonotonic group-life. 

2 Preliminaries 

Definition 2.1 (Dhaene J. [3]) Let R n be a n-dimensional Euclidean space, 
A C R n is said to be a comonotonic vectors set if for any vectors x, y € A, it 
holds that x ^ y or x ^ y . Where x ^ {^)y means that xi ^ (^)yi, i = 1, . . . ,n. 

Following 3 s , we see that the earlier concept of two comonotonic random 
variables due to Schmeidler (1986), Yarri(1987) 's work on economics. Let 
(X, Y) be a two dimensional random vector on the probability space (ft, J 7 , P). If 
for any wi,aj2 6 ^ there always holds that(X(wi) — X(u}2))(Y(u\) — Y(lj2)) ^ 0, 
then, (X, Y) is said to be comonotonic, this definition was weaken as (X(uji) — 
X{uj2)){Y(uji) — Y(u>2)) 3* 0,P.a.s. around 1997. For a n-dimensional random 
vector X = (Xi, . . . , X n ) defined on (f2, J 7 , P), if there exists a subset A C R n 
such that P(X g A) = 1, then A is said to be a support of X. A Random vector 
X = (Xi, . . . , X n ) is said to be comonotonic if its support is comonotonic. 
Lemma 2.1 (Dhaene J. et al. 3 ) A random vector X = (Xi, ■ ■ ■ ,X n ) on prob- 
ability space (JljJ 7 , P) is comonotonic if and only if the following equivalent 
conditions hold: 

(1) X has a comonotonic support; 

(2) For all values of X, x = (xi, • • ■ , x„), the joint distribution of X is 

Fx(x) = min{F Xl (a;i),Fx 2 (x2), . . .,F Xn (x„)}; 
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(3) Let U be a random variable uniformly distributed on interval (0,1). Then, 

X= d (F^(U),F^(U),...,F^(U)); 

(4) There exits random variable z and non-decreasing real valued functions fi(i = 
1 n) such that X= d {fi(z)J 2 (z),...,f n (z)). 

Following the example of pQ , [5] , [5] , the dependent multiple decrements mod- 
els for single life in life insurance have been investigated extensively. Here we 
only introduce some basic notation and main proposed methods, for other de- 
tails the readers are referred to [T],[2],[S] and the literatures therein. 

We consider a group of lives, exposed to m competing causes of withdrawal 
from the group. It is assumed that each individual may withdrawal from any 
single one of the m causes. It is assumed that , at birth, each individual is 
assigned a vector of times T%, ■ ■ ■ , T m , ^ Tj < oo, j = 1, • • • , m, representing 
individual's potential lifetime, if the individual were to withdrawal from each 
one of the m causes. Obviously, the actual lifetime span is the minimum of all 
the Ti, • • ■ , T m . Thus, it is clear that under this model the lifetimes T\, ■ ■ ■ , T m 
are unobservable, and we can only observe the min(Ti, • • ■ , T m ). In the classical 
multiple-decrement theory the random variables Xi, • • • , T m are assumed inde- 
pendent, however, decrement in many real-life actuarial applications tend to be 
dependent, and the random variables Ti, • • • ,T m are considered stochastically 
dependent and also non-defective, i.e. P(Tj < oo) = 1. The joint distribution 

F(t\, • ■ ■ , tfn) — P(Ti < ii, ■ • • , T m < t m ), 

the multivariate joint survival function 

S(ti, ■ ■ ■ , t rn ) = P(T\ > t%, ■ ■ ■ , T rn > t m ) 

which is considered absolutely continuous and where tj > 0, for j — 1, • • • ,m. 
The overall survival function of an individual aged x > is defined through 
random variable min(Ti, • • • , T m ) as 

S(i) := S{t, ■■■,t) = P(Ta >t,---,T m >t) = P(min(Ti, • • ■ , T m ) > t). 

The crude survival function (t) is defined as 

(t) = P(min(T x , . . . , T m ) > t, min^, • • • , T m ) = Tj). 
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The net survival function S'^(t) is defined as S'W(t) = P(T 3 > t). Note that 
S'^(t) is the marginal survival function, due to cause alone, associated with 
the joint multivariate survival function 

S(ti , • • • , t m ) = P(Ti > t\ , ■ ■ ■ , T m > t m j. 

Thus, we can view F(t\, • • • , t m ) = P(T\ < t±, ■ ■ ■ , T m < t m ) as a multivariate 
distribution with marginal distributions F ^'(t) = 1 — S ^'(t), j = 1, • • • ,m. 
If we know S ^'(i), we can identify and calculate the joint survival function 
S(t\, • • • , t m ) 1 and hence evaluate the overall survival function S(t, ■ ■ ■ , t) under 
some assumption of copula functions. 

The copula is one of the most useful tools for handling multivariate distri- 
butions with given univariate marginals Fx , • • • , F m . Formally, a copula C is 
a cumulative distribution function, defined on [0,l] m , with uniform marginals. 
Given a copula C, if one defines 

F(xi,...,x m ) = C(F 1 (x 1 ),F 2 (x 2 ),- ■ ■ ,F m (x m )),{xi,- ■ ■ ,x m ) G R m , 

then F is a multivariate distribution with univariate marginals Fx , ■ • • , F m . Ac- 
cording to Sklar theorem (1959) [5], Given a continuous joint distribution func- 
tion F(xx, • • • , x m ) of m-dimensional random vector X = (Xx, ■ • ■ , X m ), with 
marginal distribution functions 

Fi(xi),F 2 (x 2 ), ■ ■ ■ ,F m (x m ), (xx, ■ ■ ■ ,x m ) S R m , 

there corresponds to it a unique m-dimensional Copula function C can be con- 
structed as 

C( Ul ,---,u m ) =F(Fi 1 (u 1 ),-- .,J'- 1 (u m )),(u ll . ■-,u m ) e [0,l] m . 

Note that different multivariate distributions F may have the same copula. 
Most of the multivariate dependence structure properties of F are in the copula 
function, which is independent of the marginals and which is, in general, easier 
to handle than the original F. 

The Sklar theorem can be restated to express the multivariate survival func- 
tion S(xx, • ■ • , x n ) via an appropriate copula C called the survival copula of 
{Xx,---,X n ). Thus, 

S(xx, •••,!„) = C(Sx(xx), • • ■ , S n (x n )). 
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Where Si(xi) = 1 — Fi{xi) = 1 — Ui,i = 1, •••,n. When n = 2, we have 
C(l — 1 — ^2) = 1 — ui — i«2 + C(Tii,Ti2)- There are many copula functions 
constructed by using the given distributions, e.g. the multivariate Gaussian 
copula, Student's i-copula, and the popular Archimedean copulas constructed 
by C A {u\ 1 U2) = (4>(ui) + 0(1*2)), where <f> is a continuous, convex function 
called generator, such that 0(1) = 0,0(0) = +00, of which the frequently used 
Clayton copula, Gumbel copula, etc. (see. Nelsen (1999)0). 

Some decrement models for independent as well as comonotonic group-life 
with specified status under independent causes of decrements have been ob- 
tained by Wang ( see. Wang Z.J. and Wang D. (2009) [TO]). 

3 Model for comonotonic group-life under de- 
pendent decrement causes 

Let x\, ■ ■ ■ , x n be a group of n lives with comonotonic structure, and each mem- 
ber Xi of the group can be withdrawn with one of m different dependent causes 
of decrements T\, ■ • ■ , T m . We start with the simple case where the status of the 
lives group are determined as g(x\, • • • , x n ), e.g., which can be the joint-survival 
status (xi ■ ■ -x n ), or the last-survival status (xl ■ ■ ■ x n ), etc.. g[x\, ■ ■ ■ ,x n ) can 
be withdrawn with anyone of m different dependent causes of decrements. 

Ti(g{xi, ■ ■ ■ ,x n )), ■ ■ ■ ,T m (g(x!, ■ ■ -,x n )) 

are the latent decrement times of any status g(xi 7 ■ ■ ■ ,x n ) of the group of n 
lives. We can only observe the time 

T(g(xi, ■ ■ •,£„)) = min{Ti(3(xi, • • • ,x n )), ■ ■ ■ ,T m (g(xi, ■ ■ -,x n ))}. 

Since the dependency of the decrements causes, the latent decrement times 

Ti (^(xi, • • •,£„)), • • -,T m (g(xi,-- -,x n )) 

are dependent. Assume that Ti(g(x\, • • • , x n )), • • • , T m (g(xi, • • • , x n )) has joint 
distribution function and marginals as 

F Tl ,T 2 ,-,T m {ti,t 2 , ■ ■ ■ ,t m ) = P(T 1 (g(xi,---,x n )) ^ ti, ■ ■ ■ ,T m (g(xi, ■ ■ ■ ,x n )) ^ t m ); 
F^ih) = P{T 1 (g(xi, ■ ■ •,!„)) < ti), ■ ■ ■ ,F Tm (t m ) = P(T m (g(x 1 , ■ ■ ■ ,x n )) < t m ), 
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respectively. According to the Sklar theorem, there exists a copula function C 
such that 



FT 1 ,T 2 ,---,T m (h,t2, ■•■ , t m ) — C(Ft 1 (ti), • • • , F Tm (t m )). 

Example 1 Consider the case where g(x,y) — (xy) and g(x,y) — (xy), and 
lives (x), (y) are comonotonic and m = 2. 

(1) Let FT 1 ,T 2 (ti,t2) be the joint distribution function of T\(xy),T 2 (xy). 
Let Fr^ti) , FT 2 (t 2 ) be the distribution functions of Ti(xy), T 2 (xy) respectively. 
Then, there exists a Clayton copula function C such that 

_ JL 

F Tl .T 2 {tiM) = [(F Tl (h))- e + (F T2 (t 2 ))- e - l] \ 
where the parameter 9 > 0, and 

F Tl (ii) = max{ tl(! W, fl q^}, F T2 (t 2 ) = max{ t2 4 2 \ t2 g<, 2 )} 

by the lemma 2.1 (3). Where tqi^ ■— P(Ti(x) ^ t) and so are in the sequel. 

(2) The conditions of (1) but g(x,y) = (xy). Then, there exists Gumbel 
copula function C such that 

F Tl ,T 2 (tiM) = ex P { - [(-ln(F Tl (ti))) e + (-ln(F T2 (t 2 ))) ] ' }, 

where the parameter 9 > can be determined by the given values of the Kandcl 
statistics t$ = 1 — |, and 

F Tl (*i) - min{ tl q£\ tl q^},F T2 (t 2 ) = mm{ t2 q^ , t2 q^} 

by the lemma 2.1 (3). 

(3) For the observable decrement time T(g(x,y)), g(x,y) = (xy), we have 

t gW - P(T(xy)^t,j = 1,2) 

= P(min{Ti(a^),T 2 (;ry)}<t) 

= 1 - C(l - max{ t g«, t <?«}, 1 - max{ t 4 2 \ t q^}) , 

where C is a survival copula function which can be transformed into Gumbel 
or Clayton form. Similarly, we have 

tqS = 1 - ?(l - mm{ t qW , t 1 - min{^ 2 ), t q^}) 

for status g(x,y) = (xy). 

Note that here the decrement probabilities tqx~\t <ly~\t qx\t qy \ m fact, are 
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the net decrement probabilities corresponding to the net survival functions 
S ^(t),S ^(i), we could not obtain them by the intuitive observations. Ac- 
cording to Carriere JF. (1994) [T] and V.K. Kaishev et al. (2007)0, using a 
nonlinear system of differential equations which represent the functional rela- 
tion between the crude survival function (t) and the net survival function 
S'^it), under condition that S^\t) can be estimated based on intuitive ob- 
servations, we can solve and estimate S ^(t). Therefore, we can obtain the 
estimates of t Qx\t Qy ,t q£\t q v 2) ■ 

Consider a general case where the status of the group of n lives are not 
specified, and the group of n lives is comonotonic and each life of the group 
exposed to m dependent decrements causes. Then, we have m dependent 
comonotonic future lifetime vectors as Ti := (T^xi), • • ■ , Ti(x n )), i = 1, • • • , m. 
In fact, the conclusion that there are dependent relations between the vectors 
T\, ■ ■ ■ , T m can be easily verified by that the covariance matrix COV(Tk, Tj) ^ 
0,k ^ j (since Cov(Ti(xi),Tj(xi)) ^ 0.) Note that this dependency in gen- 
eral can be modeled by the multivariate distribution function of the vector 
(Ti, • • • ,T TO ) which has a dimension nm, we also note that each component 
vector Ti has a joint distribution F^ (I,-) which equals to a comonotonic copula 
xnm{F T .( x j(tix), ■ ■ ■ , F T ^ Xri) (t in )} by Lemma 2.1 (2), where U := • ■ ■ ,t in ). 
It is desirable to construct the multivariate distribution function F of the vector 
(Ti, • • • ,T m ) with the multivariate marginals , ■ ■ ■ ,F^ , the comonotonic 
copulas min{i ; j' 4 ( xi )(tii), • • • , F Ti ^ Xri j(ti n )}, i = 1, • • • , m. According to H.Lin et 
aZ.(1996) only the copula of independent case can be considered for modeling 
the joint multivariate distribution function with the distributions of multivariate 
marginals, i.e. if 

%,...,f m = ^\^Ti ' ' ' ' ^T m )- 

then the copula C{u\, • • • , u m ) — YiTLi u i- However, this is not our case. The 
comonotonicity makes things much easier than above case. 
Theorem 3.1 Let x\, - ■ ■ ,x n be a comonotonic group-life exposed to m de- 
pendent decrement causes T\, ■ ■ ■ ,T m . Assume that for one member xi of the 
group the future lifetime vector (Ti{x{), ■ ■ ■ ,T m (xi)) has a joint distribution 
FTi(xi),--,T m (xi) with marginals Fj^x;)' ■ • • ,Fx m / X i)- Then, their copula C de- 
rived from Sklar theorem is a copula for modelling the dependency in comono- 
tonic future lifetime vectors T±, ■ ■ ■ , T m . 
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Proof By Sklar theorem there exists a copula function C such that 

FT 1 {x l ),---,T m (x l ){hl, ■ ■ ■ ,t m l) = C{F Tl ( xl ){tu), ■ ■ ■ ,FT m {xi){tmi)), 

and F Tz(xi) (Ui) ^ mm{F T . (xi) (t a ), ■ ■ ■ , F Ti(Xn) (t ln )} = F T .(U),i = l,---,m, 
since copula C is a distribution function, thus we have 

C(F Tl{xi) (t u ), ■ ■ ■ ,F Tm{xi) {t m i)) ^ C(F Ti {h), ■■■,F T Jt m )), 

the most dependency structure properties of the vectors T\ , ■ ■ ■ , T m are in the 
later function C(Fy i (ti), ■ ■ ■ , Fj; (t m j), which can be defined as a joint distri- 
bution function 

of mn-dimension. □ 

Corollary 3.1 Assume the condition of Theorem 3.1. Then, the copula C is 
also a survival copula for the joint vector survival function . . . (t\ , ■ ■ ■ , t m ) , 
i.e. 

^Ti,---.T m ( 1 ' ' ' ' ' m ) = ^(^Ti (^l)i ' ' ' > ^T m (fm))- 

Proof By Lemma 2.1 (3), we have 

%,(*<) = P{Ti{xi) > tii, • • ■ ,Ti(x„) > ti„) - P{F T ^ xi) (U) > t a , ■ ■ ■ , ^"^(£7) > t m ) 
= P(U > F Ti{xi) {t a ),--- ,U > F T . {Xn) {t in )) =max{F T . {xi) (t il ),---,F T . {Xn) (t in )} 
> F Ti ( Xl )(tii),i = l,---,m. 

Thus, 

C(F Tl ( Xl )(tu), F Tm ( Xl )(t m i)) ^ C{Si Ti (ti), • • • , S Tm (t m )), 

the later copula can be viewed as a survival copula and we define it as the joint 
vector survival function ... f m (fi, • ■ • ,t m ). □ 

The overall survival, crude survival for the future life time vectors can be 
similarly defined. It is important to note that min{Ti, • • ■ , T m } is the ob- 
servable decrement time of whole group-life. The crude survival S^(t) := 
P(min{Ti, • • • ,T m } > t, min{Ti, • • • ,T m } = Tj). If we view the joint vector 
survival function ... <f m (£i, ■ • • ,t m ) as a mn-dimensional joint survival func- 
tion 

^T 1 (x 1 ), — ,T 1 (x n ),—,T m (x 1 ), — ,T m (x n )(tll: - ' - j tin, 1 - 1 j t m l, • ' ' , t mn ), 
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Then, the crude survival function for the case where xi is selected as a represen- 
tative of the group can be denned as (t) = P( min{Xi(xi), • • • , Ti(x n ), • • • , T m (xi), • • • , T m (x n )} > 
t, 

min{Ti(xi), • • • ,Ti(x n ), ■ ■ ■ ,T m (xi), • • • ,T m (x n )} = Tj(xi)). 
Theorem 3.2(1) If the copula C in theorem 3.1 and corollary 3.1 is differ- 
entiable with respect to jth variable and Sj;.(tj) is partially differ entiable with 
respect to tji > for all j = 1, • ■ • , m, then 

d d(S T (tA) 

±s^(t) = c j (s Ti Ct 1 ),..-,s Wm Ct m )) x ^ 3 

where Cj(ui,- ■ -,u m ) = ^-C(ui, ■ ■ ■ ,u m ). 

(2) If the copula C is n-order differ entiable with respect to its every variables and 
each survival function S^.(tj) is n-order partially differentiable with respect to 
its every variables, then 



-S^{t) = 



dtidt 2 ■ ■ ■ dt n dtjidtj 2 ■ ■ ■ dtj n 

especially when n = 2, it holds 
d 2 S^(h,t 2 ) = dS 7j d 2 C dS Tj 

dhdh ~ dt ol as* dt j2 



C(Sf i (ii), • • • , Sf m (t m )) 



8C d 2 S T ^ 



dSyp. dtjidtj 2 



(tjl,-,%n)=(tl,-,tn) 



(tj 1 ,tj 2 )—(tiM) 



(tji,tj 2 )—(ti,t 2 ) 
j = 1, • • • ,m. 

Proof (1) See the proof of theorem 6 and lemma 1 in Carriere JF.[T](1994). 
(2) Since 

S^(t) = P(mia{Ti, T m } > t, min{Ti, • • • ,T m } = Tj) 



P(T j >t j ,T k >T j (k^j)) 



OO pOO I pOO poo poo poo 



tjl 



/(ill, • • • j tin, ' " ' j imli ' ' - 7 tmn) dtki ■ • • dtfcn f <%1 ' " ' dtj n 



poo 




It! 


■/;{ 



dtjldtj 2 ■ ■ ■ dt jn 
(hi tin, ■ ■ 

dtj\dtj2 • • • dtj n 



poo 




111 


/:{ 



ST 1 (x 1 ),---,T 1 (x n ),---,T m (x 1 ),---,T m (x n ) 

idtji ■ ■ ■ dtj n 
„ . 

C(Stp i (ti), ■ ■ ■ , (t m j) 



> dtji • • ■ dtjn- 



(*fcl, — ,ijL B )=(*fi)-"i*in))Vfc - 



Thus, 



dtidt 2 ■ ■ ■ dt, 



-S^{t) = 



d (r 



dtjidtj 2 ■ ■ ■ dtjn 



-C(Sif i (ii), • • • , S7p m (t m )) 



(tjl,-,*3'n)=(tl,-,tn) 



□ 

Theorem 3.2 indicates that in the case of comonotonic group-life, the rela- 
tionship between crude survivals and net survivals is formulated by some much 
complicated n-order non-linear partial differential equations, which in most case 
is unsolvable when we have some knowledge of the crude survivals. When n = 2, 
with the conclusion (2) of theorem 3.2, If we have obtained the copula function 
C through considering some member of the group, then by this fixed copula C, 
we may use some numerical methods introduced in Carriere JF.pQ and Vladimir 
K Kaishev et aL[S] to obtain the estimates of <Sp • • • , (i m ) based on the 
observations of the crude survivals S^(t),j — 1, • • • , m. Then the estimation 
of joint survival Srj; ^ ■ ■ • ,t m ) can be carried out. 
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